FTUV/99-19, IFIC/99-20 



bep-ph/9903473 



Seasonal Dependence in the Solar Neutrino Flux 



^1- 



P. C. de Holanda 1 - 2 , C. Pena-Garay 1 , M. C. Gonzalez-Garcia 1 , and J. W. F. Valle 1 
1 Instituto de Fisica Corpuscular - C.S.I.C. 
Departamento de Fisica Tedrica, Universitat of Valencia 
461 Burjassot, Valencia, S pain 



http:/ /neutrinos, uv.es 



2 Instituto de Fisica Gleb Wataghin 
Universidade Estadual de Campinas, UNICAMP 13083-970 



Campinas, Brazil 



MSW solutions of the solar neutrino problem predict a seasonal dependence of the zenith angle 
distribution of the event rates, due to the non-zero latitude at the Super-Kamiokande site. We 
calculate this seasonal dependence and compare it with the expectations in the no-oscillation case 
as well as just-so scenario, in the light of the latest Super-Kamiokande 708-day data. The seasonal 
dependence can be sizeable in the large mixing angle MSW solution and would be correlated with 
the day-night effect. This may be used to discriminate between MSW and just-so scenarios and 
should be taken into account in refined fits of the data. 
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The difference in the v e fluxes during the day and the 
night due to the regeneration of the v e in the earth mat- 
ter - the so-called day-night effect - is one of the mile- 
stones of the MSW solutions of the solar neutrino prob- 
lem (SNP) This effect is negligible in the just-so 
picture (^] . Conversely, it is well-known that vacuum os- 
cillations lead to a seasonal effect due to the fact that the 
distance between the Earth and the Sun changes in differ- 
ent seasons of the year leading to different predicted event 
rates beyond the simple geometrical factor. Though rec- 
ognized in the early days of the MSW effect Q the sea- 
sonal effect has been neglected in most discussions of the 
MSW solution to the SNP and has even been recently 
claimed to be absent in the MSW picture |5||| . 

Recent Super-Kamiokande data after 708 days |Q ex- 
hibit an excess of the number of events during the night 
|^| . Though not yet statistically significant this provides 
some hint in favor of the possible existence of a day- night 
effect. On the other hand there is also some hint for a 
seasonal variation in these data, especially for recoil elec- 
tron energy above 11.5 MeV. While the former would be 
an indication in favour of the MSW solution, the latter 
would favour the just-so solution. 

Here we call the attention to this interesting feature 
of the MSW solution, namely that the expected MSW 
event rates do exhibit a seasonal effect due to the different 
night duration throughout the year at the experimental 
site, which leads to a seasonal-dependent v e regeneration 
effect in the Earth. Taking into account the relative po- 
sition of the Super-Kamiokande setup in each period of 



the year, we calculate the distribution of the events along 
the year both for the large mixing angle (LMA) and the 
small mixing angle (SMA) solutions to the SNP. We find 
that the effect can be as large as the one expected in the 
just-so scenario, especially in the LMA solution, where 
it amounts to ~ 10% at the best fit point for the solar 
neutrino event rates given by 0. For the SMA solution 
we find that the magnitude of the seasonal MSW effect 
is very small at the best fit point increasing as sin 2 29 
increases within the 99% CL region. We illustrate this 
behaviour in Figs, [l] and ^ and in table |[ 

Let us now describe our calculation. For simplicity, let 
us consider the two-neutrino mixing case 
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We have determined the solar neutrino survival proba- 
bility P ee in the usual way, assuming that the neutrino 
state arriving at the Earth is an incoherent mixture of 
the v\ and vi mass eigenstates. 
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where P^ n is the probability that a solar neutrino, that 
is created as v e , leaves the Sun as a mass eigenstate i>\, 
and PEarth j g probability that a neutrino which en- 
ters the Earth as V\ arrives at the detector as v e . Similar 
definitions apply to Pj,"" and Pf e arth . 

The quantity Pjj™ is given, after discarding the oscil- 
lation terms, as 
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where Plz denotes the standard Landau- Zener probabil- 
ity ||Io|| and ^ m (ro) is the mixing angle in matter at the 
neutrino production point. In our calculations of the ex- 
pected event rates we have averaged this probability with 
respect to the production point assuming the production 
point distribution given in inn. 

In order to obtain P Eart we integrate the evolution 
equation in matter assuming a step-function profile of 
the Earth matter density. In the notation of Ref. [03, 
1 - P E 



we obtain for P E e 
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where 9 is the mixing angle in vacuum and the Earth 
matter effect is included in the formulas for Z, W\ and 
W 3 , which can be found in Ref. Q . p E e arth depends on 
the amount of Earth matter travelled by the neutrino in 
its way to the detector, or, in other words, on its arrival 
direction which is usually parametrized in terms of the 
nadir angle, $, of the sun at the detector site. 

It is very important to realize that the daily range 
of variation of the nadir angle depends on the period 
of the year. As a result the quantity p E e arth is seasonal 
dependent. This will, in turn, manifest itself as a seasonal 
dependence of the expected neutrino event rates. The 
general expression of the expected signal in the presence 
of oscillations at a given time t, S osc (t), is 



S osc {t) = JdE„\{E„)x [a e {E v )P ee {E v ,t) 
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-a x (E v )(l-P ee (E v ,t))], 



where E v is the neutrino energy, A is the neutrino energy 
spectrum |l3| ] with the latest normalization pif , a e (a x ) 
is the v e (i/ x , x = /x, r) interaction cross section in the 
Standard Model |l5| , and P ee is the v e survival probabil- 
ity, which varies in time through the interval of day and 
night along the year. The expected signal in the absence 
of oscillations, S no ~ osc , can be obtained from Eq.(^) by 
substituting P ee = 1. 
The cross sections 

(Tex <W6 calculated including radia- 
tive corrections and must be corrected for energy thresh- 
old and resolution effects. In the calculation of the ex- 
pected signal it is understood that the u a -e cross sections 
a a (E) (a — e, x) have to be properly corrected to take 
into account the detector energy resolution and the anal- 
ysis window for each experiment. In Super-Kamiokande, 
the finite energy resolution implies that the measured ki- 
netic energy T of the scattered electron is distributed 
around the true kinetic energy T' according to a resolu- 
tion function Res(T, V) of the form fu|: 
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where 



s y/T'/MeV 



(6) 
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and s = 0.47 MeV for Super-Kamiokande 0n tne 

other hand, the distribution of the true kinetic energy T 1 
for an interacting neutrino of energy E v is dictated by 
the differential cross section da a {E V: T')/dT', that we 
take from Iffal. The kinematic limits are: 



< T' < T (E v ) , T (E v ) 
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For assigned values of sq, T m j n , and T max , the corrected 
cross section a a (E v ) is defined as: 
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Finally, in order to compare our results with the recent 
data from Super-Kamiokande collaboration, we must also 
include the geometrical seasonal neutrino flux variation 
due to the variation of the Sun-Earth distance [L w 
1.5 x 10 13 cm) arising from the Earth's orbit eccentric- 
ity because the neutrino fluxes in Eq. (||) are yearly aver- 
ages. In order to account for this effect we assume a 1/L 2 
dependence of the flux. Notice that Super-Kamiokande 
data are presented as ratio of observed events over the 
expected number in the Standard Solar Model where this 
expected number of events does not include the geomet- 
rical variation. Thus we must compare the experimental 
points with the predictions: 



where 
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and e = 0.0167 is the eccentricity of the Earth's orbit 
around the Sun, and T = 1 year. 

We now turn to our results. In order to study the 
behaviour of the seasonal variation we have explored the 
parameter space around the small and large mixing angle 
solutions, SMA and LMA, respectively. We find that 
depending on the values of the mass and mixing angle, 
one may get either an enhancement or a damping of the 
geometrical effect or even more complicated variations of 
the signals. 

In Fig. Q we present the expected event numbers 
in the recoil electron energy range from 11.5 MeV up 
to the maximum, plotted versus the period of the year 
for different points in the SMA solution region of the 
SNP divided by the BP98 SSM predictions in the ab- 
sence of neutrino conversions O] . We plot the expected 
behaviour for three points: the best fit point obtained 
by P| with an arbitrary S B flux, Am 2 = 5. x 10 ~ 6 eV 2 
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FIG. 1. Ratio of predicted event rate to the SSM predic- 
tion versus time of the year in Super-Kamiokande for various 
points in the SMA solution region of the SNP as labelled. We 
have normalized these three curves to the same yearly aver- 
aged event rate which corresponds to 8 B flux normalization 
0.7 for the best fit point. We also show the expectation in 
the absence of oscillations with 8 B flux normalization of 0.47 
(short dashed line) and the expected effect in case of vacuum 
oscillations for solution C in Ref. |j| (dash-dotted curve) . 
The 708 Super-Kamiokande data points are also displayed. 

and sin 2 29 — 3.5 x 10~ 3 , a point inside the 99% confi- 
dence level allowed region with Am 2 = 8. x 10 -6 eV 2 and 
sin 2 29 — 8.x 10~ 3 and a near point with Am 2 = 8. x 10~ 6 
eV 2 and sin 2 29 = 1.2 x 10~ 2 . We have normalized these 
three curves to the same yearly averaged event rate. This 
corresponds to a S B flux normalization 0.7 for the best 
fit point. For the sake of comparison we also plot the ex- 
pected behaviour in the absence of oscillations with 8 B 
flux normalization of 0.47 as well as the best fit point for 
the vacuum solution C of Ref. || . 

As seen in the figure the seasonal effect is comparable 
to the expectation in the absence of oscillation at the 
best fit point of the SMA solution and it increases as the 
mixing angle increases. In Table || we show the seasonal 
variation (in percent) defined as 



Var 



. Rn 



Rn 



+ Rn 



for the different MSW and vacuum solutions of the SNP 
where R(t) = N osc (t) / N SS m- We find that for the SMA 
solution the effect increases as one increases sin 2 29. For 
example for sin 2 29 = 0.008, still within the 99 % CL 
allowed region, it reaches 10% and for sin 2 29 = 0.012 it 
gets to be as large as 20%. Of course, since the seasonal 
effect is induced by the variation of the regeneration in 
the Earth along the year, the effect is large only in the 
parameter region where the day-night effect is not neg- 
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TABLE I. Seasonal variation (in percent) of the ratio of 
predicted event rate in various oscillation scenarios to the 
SSM prediction. 



ligible, which corresponds to larger mixing angle values. 
Note that in the SMA region the points we have chosen in 
order to illustrate the possible seasonal variation in the 
MSW picture are consistent with the measured yearly 
average day-night asymmetry. 

Now we turn to the LMA solution of the SNP where 
the effects are potentially larger. Our results for this case 
are displayed in Fig. ||. Again, we plot the expected be- 
haviour for three characteristic points: the best fit point 
obtained by || with an arbitrary 8 B flux (sin 2 29 = 0.57, 
Am 2 = 1.6 x 10" 5 eV 2 ), a point inside the 99% confi- 
dence level allowed region with Am 2 = 1. x 10 -5 eV 2 
and sin 2 29 — 0.6 and a point inside the allowed re- 
gion where the expected average day-night asymmetry 
is smaller, Am 2 = 3.2 x 10~ 5 eV 2 and sin 2 29 = 0.6. We 
have normalized these three curves to the same yearly 
averaged event rate. This corresponds to a S B flux nor- 
malization 1.45 for the best fit point. We also plot the 
expected behaviour in the absence of oscillations with 8 B 
flux normalization of 0.47 and the best fit vacuum solu- 
tion C from Ref. §. In Table | we show the variation 
(in percent) corresponding to these points. As seen in 
the table the effect at the best fit point of the LMA solu- 
tion (10 %) is comparable with the corresponding effect 
in some of the favoured vacuum oscillation solutions. In 
the LMA solution region the seasonal variation is very 
mildly dependent on the mixing angle while presents an 
oscillatory variation with Am 2 . Our results show that 
depending on the mass and mixing angle values one may 
get an enhancement or damping of the geometrical effect. 
We must bear in mind, however, than in the lower Am 2 
part of the LMA solution region, the expected yearly av- 
erage day-night asymmetry is in conflict with the existing 
data. 

Finally let us comment on the effect of an enhanced 
hep neutrino flux as suggested in p8| in order to ac- 
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FIG. 2. Ratio of predicted event rate to the SSM predic- 
tion versus time of the year in Super-Kamiokande for various 
LMA solutions of the SNP labelled in the figure. These three 
curves are normalized to the same yearly averaged event rate 
corresponding to a S B flux normalization 1.45 for the best 
fit point. We also show the expectation in the absence of 
oscillations with 8 B flux normalization of 0.47 (short dashed 
line) and the expected effect for vacuum oscillation solution 
C in Ref. Q (dash-dotted curve) together with the 708 Su- 
per-Kamiokande data points. 

count for the recent Super-Kamiokande measurements of 
the energy spectrum. We find that even with large hep 
enhancement factors of 20 or more, the expected modi- 
fications of our results near the best fit points both for 
the SMA and LMA are small. 

To summarize, we have shown that MSW solutions of 
the solar neutrino problem can lead to a sizeable sea- 
sonal dependence of the event rates in the large mixing 
angle region and this should be taken into account in 
refined fits of the data where the day-night analysis is 
also performed. The MSW seasonal effect is correlated 
with the day-night asymmetry and may potentially be 
useful in order to pinpoint the underlying mechanism in- 
volved in the explanation of the solar neutrino anomaly, 
discriminating between different solutions. For example, 
the non-observation of the day-night effect and the con- 
firmation of seasonal-dependent rates would provide an 
indication for the just-so picture. Conversely, a possible 
confirmation of a seasonal dependence accompanied by 
the day-night effect would point towards a LMA MSW- 
type solution. 
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